M = I Wess Zumino Model in ci = 3 at zero and finite temperature 
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Supersymmetric renormalization group (RG) flow equations for the efi^ective superpotential of 
the three-dimensional Wess-Zumino model are derived at zero and non-zero temperature. This 
model with fermions and bosons interacting via a Yukawa term possesses a supersymmetric ana- 
logue of the Wilson-Fisher flxed-point. At zero temperature we determine the phase-transition 
line in coupling-constant space separating the supersymmetric from the nonsupersymmetric phase. 
At flnite temperature we encounter dimensional reduction from 3 to 2 dimensions in the infrared 
regime. We determine the finite-temperature phase diagram for the restoration of the global Z2- 
symmetry and show that for temperatures above the phase-transition the pressure obeys the 
Stefan-Boltzmann law of a gas of massless bosons in 2 + 1 dimensions. 
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I. INTRODUCTION 



In this paper we investigate the three-dimensional 
Af = 1 Wess-Zumino model with general superpotential 
and explore the model beyond the realm of perturba- 
tive expansions. This quantum field theory describes 
Majorana fermions and uncharged bosons in interaction, 
with their spatial motion restricted to a two-dimensional 
layer. The selfcoupling of the bosons and the Yukawa 
coupling between fermions and bosons are such that the 
theory possesses one supersymmetry. There exists a class 
of superpotentials for which the three-dimensional mod- 
els are perturbatively renormalizable, in contrast to the 
four-dimensional models. For superpotentials of the form 
W{4>) ^ XclP"''^^ there exists both a supersymmetric and 
a nonsupersymmetric phase. In this paper we shall cal- 
culate the phase-transition curve separating the super- 
symmetric from the nonsupersymmetric phase. Besides 
supersymmetry the action is invariant under — > — 0, 
and at zero temperature the breaking of this global Z2 
symmetry is intimately linked to the breaking of super- 
symmetry. We shall see that there exists a finite phase- 
transition temperature at which Z2 symmetry is restored, 
independent of our choice for the couplings at the cutoff- 
scale. Similarly as for other two-dimensional systems, 
e.g. surface science, heterostructures or electron gases, 
the physics in two space-dimensions is rather different 
from that in three space-dimensions. 

We employ the functional renormalization group (RG) 
to calculate the phase structure at zero and finite tem- 
perature, the scaling behavior of the mass with the RG- 
scale, the wave-function renormalization, critical expo- 
nents, the effective potential and the temperature depen- 
dence of the pressure. The method has previously been 
applied to a wide range of nonpcrturbative problems such 
as critical phenomena, fermionic systems, gauge theories 
and quantum gravity, see [H-Q for reviews. A number 
of conceptual studies of supersymmetric theories has al- 
ready been performed with the functional RG. The del- 
icate point here is, of course, the construction and use 
of a manifestly supersymmetry-preserving regulator. For 



the four-dimensional Wess-Zumino model such a regu- 
lator has been presented in Recently, general 
theories of a scalar superfield including the Wess-Zumino 
model have been investigated with a Polchinski-typc RG 
equation in yielding a new approach to supersym- 
metric nonrenormalization theorems. A Wilsonian effec- 
tive action for the Wess-Zumino model by perturbati vely 
iterating the functional RG has been constructed in [l^l . 

The present study builds on our earlier results on two- 
dimensional supersymmetric field theories at zero tem- 
perature (isl . [ij] as well as on supersymmetric quantum 
mechanics, where we have constructed a manifestly su- 
persymmetric functional RG flow, see [l^. The two- 
dimensional models possess an infinite series of fixed 
points described by two-dimensional super-conformal 
theories. On the contrary, supersymmetric models in 
three dimensions possess just one fixed point, similarly 
as three-dimensional 0{N) models, see e. g. [2!. [l6l - [l8| . 

In the present work we first sketch the derivation 
of the manifestly supersymmetric flow equations in 
Sect, mil Since there exist no Majorana fermions in 
three-dimensional Euclidean spacetime we derive the flow 
equations in Minkowski spacetime and continue the result 
to imaginary time. We investigate the flow of the super- 
potential in Sect. lIIl"Sl study the fixed-point structure in 
detail, and identify the supersymmetric analogue of the 
Wilson-Fisher fixed point of three-dimensional bosonic 
0{N) models with one unstable direction. Taking into 
account a nonzero anomalous dimension in Sect. IIIIBI 
yields a scaling relation between the critical exponent of 
the unstable direction and the anomalous dimension. In 
addition we determine the zero-temperature phase dia- 
gram for spontaneous breaking of supersymmetry. 

The second part of this paper is devoted to the be- 
havior of the model at finite temperature. The fate of 
supersymmetry at finite temperature has been discussed 
extensively in the literature. For example, in previous 
works the KMS condition has been implemented directly 
in thermal superspace [l^ • In [2^ l2l| the supersymme- 
try breaking has been studied on the level of thermal 
Green functions. The breaking of supersymmetry by fi- 
nite temperature corrections, for example the one-loop 
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corrections to fermionic and bosonic masses, has been 



determined in the real-time formulation in [22[. The in- 
evitable breaking of supersymmetry at finite temperature 
has sometimes been called spontaneous collapse of super- 
symmetry [23| . 

In Sect. IIVI we derive the RG flow equations at finite 
temperature. In addition to the momentum integrals we 
are confronted with sums over Matsubara frequencies. 
For the three-dimensional Wess-Zumino model and for a 
particular regulator the thermal sums can be calculated 
analytically. Related sums have been discussed in earlier 
works on finite-temperature renormalization group flow 
equations, for example in (23 - [27l [29l - l32j . We observe 
that the Wess-Zumino model in three dimensions at finite 
temperature in the symmetric phase behaves similarly 
to a gas of masslcss bosons. In particular we show in 
Sect. IIV Al that it obeys the Stefan-Boltzmann law in 
three dimensions. For high temperatures the fermions 
do not contribute to the flow equations since they do not 
have a thermal zero-mode. On the other hand we observe 
dimensional reduction in the bosonic part of the model 
due to the presence of a thermal zero-mode. We show in 
Sect. IIV Bl how this is manifested in our RG framework. 
In a similar way dimensional reduction has been observed 
in 0(iV)-models at finite temperature in [H,!!!]. Finally 
we compute the phase diagram for the restoration of the 
global Z2 symmetry at finite temperature in Sect. IIV CI 



II. THE = 1 WESS-ZUMINO MODEL IN 
THREE DIMENSIONS AT T = 

There are many works on the supersymmetric Wess- 
Zumino models in both four and two space-time di- 
mensions. Actually the two-dimensional model with 
N = 2 supersymmetries is just the toroidal compactifi- 
cation of the four-dimensional A/" = 1 model. The three- 
dimensional model with M — 1 supersymmetry, on the 
other hand, cannot be obtained by dimensional reduction 
of a local field theory in four dimensions. Thus it may be 
useful to recall the construction of the three-dimensional 
model starting from the real superfield 



The anticommutator of two supercharges yields 
{Qa,Q^} = 2{'^'^)a^ d^. The supercovariant derivatives 
are 

Pi f) 

V^— + i{-1^a)df,, and V ^ i{a-i^)d^,. [A) 

Up to a sign they obey the same anticommutation rela- 
tion as the supercharges 

{2?„,P^} = -2(7)/a^. (5) 

As kinetic term we choose the D term of V^V^ = 
2Q:a>Ckin -I- . . . which reads 

^kin = ^d^cl^d^^cj, - ^^^V + ^F'- (6) 

The interaction term is the D term of 2VF($) = aaCint + 
. . . and contains a Yukawa term, 



Ant = FW'{c^) - ^W"{(l))^^. 



(7) 
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<I'(x, a) = (j^ix) + aip{x) H — aaF{x) 



(1) 



The complete off-shell Lagrange density CoS = 'Ckin+'Cint 
takes then the simple form 

Cos = ^d^cl>d^(l>-'-4>0ij + ^F^ + FW'{(f>)~^W"{^)i;i;. 

(8) 

Eliminating the auxiliary field via its equation of motion 
F ~ —W'{(f>), we end up with the on-shell density 

Con = \d^<j,d^'cp - - \W'\^) - iw^"(0)^^. 

(9) 

From this expression we read off that W''^{(t)) acts as 
a self-interaction potential for the scalar fields. For a 
polynomial superpotential W{(j)) in which the power of 
the leading term is even, W^(0) = C(/)^" -|- 0(0^"), we do 
not observe supersymmetry breaking in our present non- 
perturbative renormalization group studjQ. On the other 
hand spontaneous supersymmetry breaking is definitely 
possible for a superpotential in which the power of the 
leading term is odd. In the explicit calculations we shall 
use a Majorana representation for the 7-matrices, 7*^ = 
(72, 1^ ^ i<73 and 7^ = tci. 



with real (pseudo)scalar fields 4>, F and Majorana spinor- 
ficld ip. The supersymmetry variations are generated by 
the supercharge 



^ da 



(2) 



We use the metric tensor (77^;^) = diag(l, —1 — 1) to lower 
Lorentz indices. With the aid of the symmetry relations 
for Majorana spinors ipx = xi^i V'T^'X = ~X7^V' £ind the 
particular Fierz identity aa = —aa 1/2 the transforma- 
tion laws for the component fields follow from Eq. ([2]): 



P^, 6'ip = {F + i 



SF = il3$ij . (3) 



III. FLOW EQUATION AT ZERO 
TEMPERATURE 

To find a manifestly supersymmetric flow equation in 
the off-shell formulation we extend our earlier results on 



^ In a two-loop calculation a ground state with broken super- 
symmetry has been found in Ref. [36| . Since we neglect higher 
_F-tcrms in our non-perturbative study it is not possible to check 
whether the findings of this perturbative analysis of the Wess- 
Zumino model hold when higher-order corrections are taken into 
account. 
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the one- and two-dimensional Wess-Zumino models [l3|, 
[Tsj to three dimensions. Since there are no Majorana 
fermions in three-dimensional Euclidean space we begin 
with a Minkowski spacetime formulation of the Wetterich 
equation [btI. Issj: 



dtTk 



-Str 



dtRk 



t = In/c^ 



(10) 



where the scale-dependent effective action interpo- 
lates between the microscopic (classical) action Tk=K = S 
and the full quantum effective action Tk=o = T- The sec- 
ond functional derivative in Eq. (fTO|) is defined as 



(r?'). 



a ^(5*6 



where the indices a, h denote the field components, in- 
ternal and Lorentz indices, as well as space-time or mo- 
mentum coordinates, i.e., ^I*"^ ~ {(j), F, ip, -0) is a vector of 
fields, not to be confused with a supcrficld. The cutoff 
function Rk provides an infrared (IR) cut-off for all fields 
and specifics the Wilsonian momentum-shell integrations 
such that the flow of Tk is dominated by modes . 
For a derivatior0 of the RG flow equation in Minkowski 
space-time (HU]) we refer to App. [X] 

To construct a supersymmetric flow we use as regulator 
an invariant D term in superspacc. Since such a term 
should be quadratic in the ficldso, see App.[^for details, 
it is the D term of a superfield with K being a 

function of VV. Using the anticommutation relation ([5]), 
powers of 'D'D can always be decomposed into 



2n 



i-ny 



(11) 



such that an invariant and quadratic regulator action is 
the superspacc integral of 



^2 



-□)-U. 



(12) 



Expressed in component fields, we find 



F 



- I ^PRU. (13) 



In the following we neglect an additional term to the flow equa- 
tion arising from the normalization of the Gaufiian measure in 
the partition function. Including such an additional term yields 
a field-independent constant to VF' (</>). We stress that only non- 
universal quantities such as the critical temperature for the 1,2 
phase transition are affected by such a constant. However, our 
analysis of the critical dynamics at the phase boundary at zero 
and finite temperature is not affected. 

A Regulator term quadratic in the fiuctuation fields ensures that 
we eventually obtain a non-perturbativc RG equation with one- 
loop structure. 



In momentum space, zi9^ is replaced by and the oper- 
ators take the explicit forms 



Rk 







Rl 



(14) 



with 



Rf = 



and R\ = -ri - ., (15) 



where ri = ri(p^/fc^) and r2 = 7'2(p^/fc^)- Note that 
the requirement that the RG flow preserves supersym- 
metry relates the regulators R^ and R^ in the bosonic 
and fermionic subsectors. 



A. Local Potential Approximation 

We employ the following ansatz for the supersymmetric 
effective action for our study of the three-dimensional 
Wcss-Zumino model: 



1 



-F' 



+FWl[4')-7^Wl!{cj))iPi^ 



(16) 



In the following we work in the so-called local potential 
approximation (LPA) where the expectation values of the 
fields are taken to be constant over the entire volume. As 
it has been found in studies of scalar 0{N) models and 
low-energy QCD models, the LPA captures already most 
qualitative and quantitative features associated with crit- 
ical dynamics at zero and finite temperature provided the 
anomalous dimensions are small, see e. g. Refs [l, HgI - 
[isl . m, lU . For the time being we restrict our study to 
LPA. In Sect. IIIIBI we shall then discuss the running of 
the wave-function renornialization. 

In order to obtain a flow equation for the superpoten- 
tial, we project Eq. pH)) onto the terms linear in the 
auxiliary field F and integrate the resulting with re- 
spect to (f). Performing a Wick rotation of the zeroth 
component of the momentum, i.e. pg^ — > ip^ , we find 
the flow equation 



dkWk{cl>) 



d^p 9fcri(l + r2)-9fcr2(W;(0)-ri) 
(27r)3 p2(^2 + l)2 + (VK^'((/))+ri)2 • 

(17) 

In the following we choose the simple regulator functions 



ri =0, 



r2 



k 

¥\ 



i]e{e^p^), 



(18) 



for which the momentum integration in ([T7)) can be per- 
formed analytically. In the present work we do not aim 
at a study the regulator dependence. However, the reg- 
ulator dependence of functional RG flows, in particular 
with respect to critical phenomena, has been investigated 
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FIG. 1. RG flow of Wk{(l)) with the initial conditions Aa/A = 
1, a\/K^-^ = 0.02. Note that t is related to via t = Infc^ 
see also Eq. (fTU)) . 



in great detail and it has been shown that optimized reg- 
ulator functions minimizing the trucational error can be 
constructed, see e. g. Refs. 0, HI [13, [H, [13^ 
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Contrary to the model in two dimensions [13[ , the reg- 
ulator function regularizes the flow even if we allow 
for running wave function renormalizations. For the su- 
perpotential Wk we then obtain the simple flow equation 



dkWu 



(19) 



As we are interested in the bosonic potential = 
PF'^((/))/2 we will consider mostly the flow equation for 
which reads 



87r2 + Wl!{(PYy 



(20) 



Fig. [T] shows the flow of W^{4)) for a quadratic super- 



potential at the cutoff scale, 
with initial conditions Aa/A 



Q, and 



1, ai/A°-5 = 0.02. With 



these initial conditions the RG flow starts in the regime 
with broken Z2 symmetry and for fc — > ends up in the 
regime with restored Z2 symmetry. We observe that the 
potential Vk becomes flat at the origin as k is lowered to 
the infrared. In addition the function W^^q is regular 
for all values of the field, in contrast to the situation in 
two dimensions. 

In order to study the fixed-point structure we introduce 
dimensionlcss quantities 



and its fixed points are characterized by dtw^ = 0. The 
fiow equations in two and three dimensions have almost 
identical forms. In three dimensions, however, there ap- 
pears the additional term oc ipw'j^^ip), since the field (f> 
itself is a dimensionful quantity. 

We observe a further peculiarity of supersymmetric 
Wcss-Zumino models: Only the second derivative of the 
superpotential enters the fixed-point equation following 
from dtw* = 0, see (T^j. It follows that the couplings of 
the terms and (j)^ do not enter the fixed-point equation 
but evolve independently. As we shall see below, this has 
some interesting consequences which distinguish the su- 
persymmetric Wess-Zumino model from purely bosonic 
theories, for ex amp le 0{N) models in three dimensions, 
see e.g. 

For our fixed-point analysis, we study the first deriva- 
tive of Eq. dUl), 



Wi. 



47r2(l 



87r2(l 



„l!2\ 



(23) 

where the prime denotes the derivative with respect to 
the dimensionlcss field ip. 



Polynomial approximation 



First we solve Eq. I\23\i in the polynomial approxima- 
tion with a Z2 symmetric w'j^ at the cutoff scale. The 
RG flow is such that a Z2 symmetric w'^ will remain Z2 
symmetric during the flow. Thus a polynomial approxi- 
mation for ui'a is of the form 



„2i 



(24) 



where A, and &2i denote the scale-dependent couplings. 
Recall that for an even function w'{lp) supersymmetry 
may be broken. We find the following infinite tower of 
differential equations: 



dtXit) = - 



3A(0^ 354(t) 



7r2 27r2A(t) 

3b4(t)-6A(t)3+7r2A(t) 
27r2 



- 1 



1 

4^ 



(25) 



dMt) = 



12064(i)A(t)2 + 2TT%i{t) - lbhfi{t) - %Q\{tf 
4^^2 



ip = k ^(j), wk{ip) = k ^l^fc(0), 
w'^{^) = k-iwl{^), .... 
The dimensionlcss fiow equation then reads 



8^2(1 +y;^'2) 



(21) 



(22) 



Note that due to supersymmetry the lowest order cou- 
pling o? does not enter the fiow equations of the couplings 
A, 64, 65, 

In our fixed-point analysis we find a Gaui3ian fixed 
point with all coupling constants equal to zero and, due 
to the Z2 symmetry, a pair of fixed-points whose cou- 
plings converge rapidly for larger truncations as shown 
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6.286 
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1.595 


2.873 


7.150 


13.41 






10 


1.595 


2.873 


7.155 


13.48 


1.212 




12 


1.595 


2.870 


7.118 


12.90 


-8.895 


-183.3 



TABLE I. Wilson-Fisher fixed point as obtained from tlie 
polynomial approximation to w'{tp). Note tfiat tfie discrep- 
ancy between the fixed-point values of &io for the truncation 
with 2n = 10 and 2n = 12 is a truncation effect. Due to our 
polynomial expansion of u)^ (</?), see Eq. (|24p . we obtain a hi- 
erachy of flow equations for the couplings b2i in which the RG 
flow of the coupling b2n (for given truncation order n) suffers 
the most by the finite truncation order n. Therefore we ex- 
pect that the fixed-point values for the higher-order couplings 
converge when n is increased, as it is indeed the case for A, 
h2, 64, be. 




FIG. 2. Ratio i)2i/i)2i+2 of the Taylor expansion coefficients 
of the potential as a function of i. 
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critical exponents 
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-0.799 


-5.92 


-20.9 
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-0.767 


-4.83 


-14.4 


-38.2 






10 


-0.757 


-4.35 


-11.5 


-26.9 -60.8 






12 


-0.756 


-4.16 


-9.94 


-21.4 -43.8 -89.0 






14 


-0.756 


-4.10 


-9.13 


-18.3 -35.1 -65.4 


-123 




16 


-0.756 


-4.08 


-8.72 


-16.4 -29.9 -52.9 


-91.9 


-163 


18 


-0.756 


-4.08 


-8.54 


-15.2 -26.4 -45.0 


-75.0 


-124 



-209 



TABLE IL Critical exponents for the Wilson-Fisher fixed 
point for different truncations. 



that for the corresponding series with dimensionful field 
and couplings 622 = • kr'^'^^'^ the radius of convergence 
shrinks with decreasing scale k, 



lim 



^21+2 



k lim 

i->cx) 021+2 



k ■ r . 



(28) 



Thus, the radius of convergence of the power series ex- 
pansion of W^(c/>) tends to zero for fc — 0. 



2. Partial differential equation 



in Tabular HI From the stability matrix, 

d{dth) 



dbi 



(26) 



we read off that the non-Gau6ian fixed points are IR sta- 
ble. Here, we have set 60 and 62 = A. These IR 
stable fixed points are to be considered as supersymmet- 
ric equivalent of the Wilson-Fisher fixed point. 

Let us now discuss the critical exponents which are the 
negative eigenvalues of the stability matrix at the fixed 
point. The coupling a^, which does not feed back into 
the equations for the higher-orders couplings, defines an 
IR unstable direction with a critical exponent l/v^= |. 
The critical exponents of the IR stable directions of the 
Wilson-Fisher fixed point are given in Tabular |lll Actu- 
ally we observe a better convergence of the lowest critical 
exponents than in two dimensions [isj . 

From Tabular U] we estimate that the radius of conver- 
gence of the Taylor series ((24|) , given by 



= lim 



^22+2 



(27) 



is finite. This can be seen more clearly from Fig. [51 where 
we plotted the ratios 621/^21+2 as functions of i. Note, 



Let us now turn to the solution of the partial differ- 
ential equation ([2^ . We have seen in Eq. ([55]) that the 
coupling associated with the IR unstable direction does 
not feed back into the fixed-point equation. It is therefore 
sufficient to consider the second derivative of Eq. ((22)) . To 
simplify the notation we introduce u{if) 
fixed-point equation for u reads 



<(<^). The 



2u- 



- 1 



/^ + 47r" ^"/'^^ {2u~^u'). (29) 
— 1 



It is straightforward to see that the above equation has an 
asymptotic solution Ui^{(p) ~ ip^ . Employing a standard 
Runge-Kutta solver for ordinary differential equations we 
find one regular odd solution for the starting condition 
u(0) = and u'(0) = 2A « ±2 • 1.59508 = ±3.19016. 
For field amplitudes ip < 0.245 the regular solution is 
bounded by it < 1 and this inner part of the solution 
corresponds to the IR stable fixed-point solution found 
in the polynomial approximation discussed above. 
Since u = w" we find for large fields 

w"{ip ±00) ~ ± =^ Wi'{(t) ±00) ~ ± (/)^ 

In other words the outer part of the regular fixed-point 
solution connects smoothly to the inner part. Thus we 
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have found a solution corresponding to a bosonic poten- 
tial V ~ (k'^^^w'^lip))^ which behaves as ^ cj/' for large 
<j>. This is the supersymmetric analogue of the Wilson- 
Fisher fixed point of three-dimensional 0{N) theories, 
see e.g. @, HI]- 



2n 


4 


6 8 10 12 


14 


V* 


0.187711 


0.188258 0.18802 0.187996 0.188001 


0.188003 



TABLE III. Dependence of the fixed-point value of the 
anomalous dimension -q on the truncation order n of the ex- 
pansion of w' in powers of the field tp, see also Eq. (|24p . 



B. Next-to- leading order 



For the next-to-leading order approximation we em- 
ploy the following ansatz 



(30) 



with Z'^ being a scale-dependent wave- function renormal- 
ization. We neglect a possible momentum and (j) depen- 
dence of Zk ■ This approximation corresponds to an inclu- 
sion of the next-to-leading order correctior0 in a system- 
atic expansion of the effective action in powers of fields 
and derivatives. As we discuss below, we find that the 
anomalous dimension i] remains small compared to one 
within this approximation, see also Tabular lllll Thus, we 
expect that higher-order corrections such as ^ {(l)d^(j))^ 
does not significantly affect our results for the zero- and 
finite-temperature phase diagran{f|. 

Projecting on the part linear in the auxiliary field and 
integrating with respect to (j) yields the superpotential. 
On the other hand projecting on the terms quadratic in 
the auxiliary fields yields the flow equation for the wave 
function renormalization. Employing the regulator func- 
tions ([T5| we find the following coupled set of differential 
equations: 



dkZl 



247r2 
47r2 



iloikdkZl 



2Zi)x 



{k^zt+w'^im 



(31) 



(32) 



the dimensionless flow equations read 

I,. . (3 - r/)rD" 



dtXc + 2tt) = -(1 + T/)xrD 
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2' ' 247r2(l 
(2-?;)(l-tr"2)ir"'2 



47r2(l + ro"2)3 



(33) 
(34) 



In order to study the fixed-point structure it is conve- 
nient to consider the anomalous dimension 77 as a free 
parameteiH, see e. g. Ref. [1^. In complete analogy to 
the two-dimensional Wess-Zumino model we find lines of 
fixed points corresponding to potentials with no nodes 
(outermost line), one node and so on, see Fig. [3] (left 
panel). In fact, we encounter exactly the same picture 
as in two dimensions, apart from a shift of the graph 
to lower 77 values. Concerning the number of fixed points 
the situation is completely different as in two dimensions: 
Because of the shift of rj we find only one pair of fixed 
points for 77 = 0, and not an infinite number of pairs. 
Such a dependence on the dimensionality has also been 
observed for 0{N) models [s^. 

Note that if we had actually used the same regulator 
as in namely 



and r2 



P 



b=0 



then the flow equation for the superpotential in two di- 
mensions would turn into the flow equation in three di- 
mensions under the transformation 7/ — > 77 — 1. This cor- 
respondence explains the similarities of Fig. [3] (left panel) 
in two and three dimensions. 

As in two dimensions, we can deduce a superscaling 
relation from the RG flow equation of which relates 
the critical exponent and the anomalous dimen- 

sion O: 



Introducing the anomalous dimension 77 
the dimensionless quantities 

x = Zuk--^^, m(x) = fc-2w^fe(0). 



-dt In Zi and 



1 



3-77 



(35) 



The truncation dependence of the fixed-point value of the 
anomalous dimension is shown in Tabular [ 



* LPA corresponds to the leading-order approximation. 

^ Note that the same reasoning has been found to hold in studies 

of the critical dynamics in 0(N) models and low-energy QCD 

models, see e. g. Refs. Q, [H, Isl, [H, [H 



Note that 77 = is a consistent solution of the fixed-point equa- 
tions. In this respect the model in three dimensions is substan- 
tially different from the model in two dimensions. There, r? = 
is not a consistent solution of the fixed-point equations. 



7 




FIG. 3. Left panel: Lines of fixed points in the 77-7 plane (solid curves) and the anomalous dimension as a fimction of 7 = 2A 
as obtained from Eq. p4|l (dotted curve). Right panel: Phase diagram in the plane spanned by the dimensionless couplings 
specified at the cutoff scale A as obtained from truncations with n = 1 {(fP'), n = 2 ((?i*), and n — i (0®) in Eq. H24p . 



results from the flow equation 
f(x)=0.23*x-0.82 



-9 -8.5 -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 

ln(k) 

FIG. 4. Logarithm of the boson mass as a function of the RG- 
scale k. A linear fit to the data points yields m(fc) ^ kP'^'^ for 
fc « A. 



C. Zero-temperature phase diagram and scaling of 
the mass 



The breaking of supersymmctry is driven by the un- 
stable direction a^, similarly as in two dimensions [isj . 
In the plane spanned by the values of the dimensionless 
couplings A and given at the cutoff-scale A we find 
a transition line for supersymmetry breaking, see Fig. [3] 
(right panel). 

From the effective potential V = limfc_j.o we 
read off the mass of the boson, 



the mass is given by 

"^broko„ = W/^^o(OX!.o(0), <_.o(0)=0, (37) 



whereas for unbroken supersymmetry it is given by 



unbroken 



(38) 



In the broken phase a polynomial expansion of the su- 
perpotential is justifieqj and we find 



"^Lko„(fc) = w^fe(oX'(o) 



From the scaling of the couplings for fc — > 0, 
A, ^ const, and ^ k~^ , 



we read off that the mass scales as 



m{k) ^ for fc < A . 



(39) 



(40) 



(41) 



This scaling behavior is demonstrated in Fig.|3]for a trun- 
cation with ri = 4, i. e. an expansion of w' up to order 
tf^ ~ (/)^, see also Eq. ((24)) . Due to our analysis of the 
convergence of the fixed-point values of the couplings A 
and &2j as well as the anomalous dimension 77*, see Tab- 
ular U [n] and mil we expect that the truncation order 
71 = 4 is already sufficient to capture qualitatively and 
quantitatively most of the features of the zero- and finite- 
temperature phase diagram, see also Sect. IIVBI From a 
linear fit to the double-logarithmic plot of m(fc) we find 
m(k) ^ fc"-^^ which is indeed reasonably close to the 
expected scaling given in Eq. (HI]). The reason for this 
behavior, which is very different from the one found in 



^'0)+W^n0(</>0): 



(36) 



where the field 4>o minimizes V . In the case of broken 
supersymmetry (and unbroken Z2) we have (pQ — and 



Note that for unbroken supersymmetry a polynomial expansion 
is bound to fail and one needs to solve the partial differential 
equation in order to determine the boson mass. 
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0{N) models, is that the unstable direction does not 
feed back into the fixed-point equatiorQ. Independent of 
the value of the coupling at the UV (ultraviolet) cutoff 
scale A, the second derivative of the superpotential flows 
always into its IR fixed point corresponding to a confor- 
mally invariant theory. In 0{N) models, on the other 
hand, the unstable direction in the broken regime feeds 
back into the fixed-point equations for the couplings. Ap- 
proaching the IR fixed point of 0{N) models requires 
therefore fine tuning of the UV parameters [3 ■ 



IV. FLOW EQUATIONS AT FINITE 
TEMPERATURE 



In this section we study finite temperature effects in 
the three-dimensional Wcss-Zumino model. To this end 
we restrict ourselves to the LPA which we expect to 
provide already a quantitative insight into the finite- 
temperature phase structure as shown in Ref. [32| . 

The finite-temperature flow equation in LPA can be 
obtained straightforwardly from the zero-temperature 
equations (|20p by replacing the momentum integration 
in time-like direction by a summation over Matsubara 
frequencies: 



to 



f(x)=1.49*x-0.66 
results from the flow equation 



ln(T) 



FIG. 5. Double logarithmic plot of the temperature depen- 
dent minimum of the superpotential versus the temperature. 



'superpotential' describing the Yukawa-type interaction 
between fermions and bosons originates in the supersym- 
metry breaking induced by the different thermal bound- 
ary conditions for the bosonic and fcrmionic flelds. 




dpo 
2tt 



(42) 



with frequencies w„ = 2'KnT for bosonic fields and Vn = 
(2n + 1)tiT for fermionic fields. We refer to Appendix iBl 
for a detailed derivation of the finite-temperature flow 
equations. Here we simply note that we can perform 
the Matsubara sums explicitly for the regulator functions 
([T5| . The flow equations read 



,hos _ ::_y^w 



(p 

(2SB + 1)' 



(43) 



k ^ ' k 

form k^ (fc^ - W'^^)Wi!' 

87r2 (fc2 + iy^'2)2 



2(2sB + l 



1-1- 



2sF7rrN2 



(44) 



where the temperature- dependent floor functions sb and 
sp are given by 



SB 



2'kT 



and sp 



27rr 



(45) 



The differences in the flow equations for the 'superpo- 
tential' describing the self-interaction of bosons and the 



Note that this is only true for finite TV. In the large TV limit, the 
running of the higher-order couplings is also independent of the 
running of the vacuum expectation value of the field Il6l| which 
corresponds to a? in our study of the M = 1 Wess-Zumino model. 



A. Pressure 

In the previous section we have shown that the boson 
mass tends to zero for A: — )■ in the phase with bro- 
ken supersymmctry (restored Z2 symmetry). Thus we 
expect that the thermodynamic properties of the three- 
dimensional model in the phase with restored Z2 symme- 
try are similar to that of a gas of massless bosons. The 
pressure of such a gas in 2 + 1 dimensions is given by 



03) r3 
27r 



(46) 



The (normalized) pressure Ap for a given temperature is 
determined by the negated difference of the interaction 
potential for a given temperature evaluated at the min- 
imum and its corresponding zero-temperature value. In 
our model this corresponds to the height of the barrier in 
the double-well potential {W'Y for a theory with broken 
Z2 symmetry or the minimum of {W'Y with unbroken Z2 
symmetry. For small temperatures we therefore expect 



A(Aa2) =(Aa2)T=o - {\a')T^Q 



C(3) 
27r 



~ cxp(-0.83)T'^ 



(47) 



for the temperature dependence of the minimum of 
VF'(^). In Fig. Owe show the double-logarithmic plot of 
the minimum of the superpotential versus temperature as 
obtained from a truncation of the potential with n — A, 
I. e. an expansion of w' up to order ip^ ~ (/)^, see also 
Eq. ([M)) . The linear fit to the double- logarithmic plot 
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yields a power law Ai^Xa^) = exp(— 0.68) • T^ '^^ whieh is 
compatible with the values for a gas of massless bosons 
in 2 + 1 dimensions, as given in Eq. (jTZj) . We expect that 
deviations from the ideal gas limit are present for two 
reasons: (i) The boson mass tends to zero for fc — )- only 
but remains finite for any k > 0, see Eq. (j4l]) and (ii) 
boson self-interactions lead to deviations from the ideal 
bose-gas limit. 



B. RG flows at finite temperature and dimensional 
reduction 

For high temperatures T ^ k both floor functions in 
(|45p vanish and the flow equations simplify considerably: 



djrr/ bos 
kWL = - — ^ 



8^2 ''1^ 



112 



2ttT 



Vt/f)2 



dkWl 



I fcrm 



0. 



(48) 
(49) 



Following Ref. [33[ we rescale the quantities in the 
bosonic flow equation according to 

- VT0, Wk (0) = TWu (0) . (50) 

For Eq. pS]) this leads exactly to the zero-temperature 
flow equation of the two-dimensional model. At this 
point, however, we should stress that the theory which we 
obtain in the limit T fk 3> 1 is not the supersymmetric 
M = \ Wess-Zumino model in two dimensions since the 
fermions have dropped out of the theory due to the ab- 
sence of thermal zero-modes. Therefore supersymmetry 
is necessarily broken at finite temperature. The fact that 
we still obtain the same functional form for the bosonic 
flow equation can be understood in terms of the role of 
the auxiliary fleld: In order to obtain the bosonic flow 
equation, we have to project on the terms coupling to 
the auxiliary field. Since there is no coupling between 
the auxiliary fleld and the fermionic part of the theory, 
the fermions do not contribute to the bosonic flow equa- 
tion. 

Because of the rescalings of the flelds and the poten- 
tial according to Eq. ([50]) we expect that the couplings 
exhibit the following behavior for T/fc ^ 1: 



T 



1/2 



ZD I 2 \ 2D( 2 \ I ^ 

(a )t = (a )t=o I 

/rp\ -1/2 

3D x 2L> \ M 

At — AT=o -r 



(51) 



^^(&2.)t='^(62.)t=o(^ 



1/2-i 



where ^^{a'^)T=o, ^^At=o and ^^{b2i)T=o denote the 
fixed-point values of the couplings of the two-dimensional 
theory. Indeed we observe such a running of the couplings 
for T/k ^ 1 in our numerical studies. 

From Eq. ([5^ we deduce the radius of convergence of 
our expansion of the potential in powers of the flelds. 



Recalling the relation between the dimensionful and di- 
mensionless couplings, b2i = 621 • k^^^~'^, we find 



3D 



3D 



r = lim 



(b2^h 



3^(62 



'2i+2jT 

^^(MT=o-(r/fc)^/^-'-fc^/^-' 

■'^^^''{b2^+2)T=0■{T/k)-y^-^■k^/^-^ ^ ' 

Thus we find a finite radius of convergence for the poly- 
nomial expansion of the superpotential at finite temper- 
ature which is a consequence of the finite radius of con- 
vergence found for the underlying two-dimensional model 
at zero temperature [l3| . However, the radius of conver- 
gence '^^r tends to zero for T 0, in full agreement with 
our results in Sect. IIII Al 



C. Phase diagram at finite temperature 

In this subsection we discuss the phase diagram of the 
Af = 1 Wess-Zumino model at finite temperature. 

Whether supersymmetry is broken or not at vanish- 
ing temperature depends on our choice for the couplings 
at the cutoff scale as we have discussed above. At fi- 
nite temperature we have necessarily soft supersymme- 
try breaking due to the different boundary conditions for 
bosons and fermions in Euclidean time direction. Besides 
supersymmetry, however, our theory is invariant under 
(f) ^ —(f). At vanishing temperature the breakdown of 
this Z2 symmetry is intimately linked to the question 
whether supersymmetry is broken or not. Recall that at 
T = broken K2-symmetry of the ground state implies a 
supersymmetric ground state and that restored Z2 sym- 
metry of the ground state implies broken supersymmetry. 
Even though supersymmetry is necessarily broken at fi- 
nite temperature, we shall see in the following that Z2 
symmetry of the ground state can be either broken or re- 
stored depending on the actual value of the temperature. 
Due to the relation between supersymmetry and Z2 sym- 
metry at vanishing temperature, we consider the strength 
of Z2 breaking as a measure for supersymmetry breaking 
at finite temperature. Thus we distinguish between the 
case of soft supersymmetry breaking due to finite tem- 
perature but broken Z2 symmetry of the ground state 
and the case with restored Z2 symmetry of the ground 
state at finite temperature. 

For our numerical study of the finite-temperature 
phase diagram we employ a 0*-truncation of the poten- 
tial^, i. e. we use n = 4 in Eq. (p4|) . 

In the left panel of Fig. [51 we show the Z2 phase- 
boundary in the space spanned by temperature T/A and 



^ The quality of such an approximation has been studied exten- 
sively for zero and finite temperatures llfiHlSl , |40|| . It indeed 
proves to be quantitatively useful in case one is not interested in 
a high-accuracy determination of critical exponents [TtI . [TsII . 
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FIG. 6. Finite-temperature phase diagram of the TV = 1 Wess-Zumino model. Left panel: TL^ phase boundary in the space 
spanned by temperature T/A and the value of the couplings (AA,a^) at vanishing temperature. Right panel: Slice of the Z2 
phase-boundary for fixed Aa = 0.8. 



the value of the couplings (Aa, a\) at the UV cutoff scale 
at vanishing temperature^. The lower end of the phase 
boundary on the T = plane in Fig. [HI corresponds to 
the phase-transition line which separates the phase with 
broken supersymmetry from the one associated with a su- 
pcrsymmetric ground state, see also Fig. [3] (right panel). 
Choosing couplings {X\,a'jJ at T = associated with a 
supersynimctric ground state (and broken S2 symmetry) 
we find always a second-order phase-transition tempera- 
ture at which the system enters the phase with restored 
symmetry. As discussed in Sect. IIV Al the thermo- 
dynamic properties of this Z2-symmetric phase are sim- 
ilar to the ones of a gas of massless bosons. We ex- 
pect that this finite-temperature phase transition falls 
into the Ising universality class with critical exponents 
determined by the underlying 2d theory, namely by 
the exponents of Onsager's solution of the corresponding 
lattice spin model. However, a quantitative analysis of 
the critical behavior would require the inclusion of the 
anomalous dimension in our studies. This is deferred to 
future work. 

In the right panel of Fig. [HI we show a slice of the 
Z2 phase boundary for fixed Aa — 0.8. We find that 
the phase-transition temperature increases with increas- 
ing (a\)T=o- Since (a^^Q)T=o increases with increas- 
ing (a^)T=o, the phase-transition temperature increases 
with increasing (ai_>.o)T=o, i-e. in terms of the renor- 
malized quantitjo. Our observation of an increasing 
phase-transition temperature with increasing (a^^Q)r=o 



As initial conditions for the finite-temperature RG flows we use 
the zero- temperature values of the couplings ( A a , ) at the cut- 
off scale. Therefore we have to restrict our study to temperatures 
significantly smaller than the UV cutoff in order to ensure that it 
is justified to use these zero-temperature values of the couplings 
at the UV starting point of the flow. 

Recall that A flows into its fixed-point value A* for T = in- 



is in accordance with our expectations from a scalar 
0(1) ~ K2 model since (A*a^_^Q)T=o sets the scale at 
T — and therefore plays a similar role as a finite expec- 
tation value of the fields in 0{N) models, see e. g. [34| . 



V. CONCLUSIONS 

In this work we have employed the functional RG for 
a study of the three dimensional A/" = 1 Wess-Zumino 
model at zero and at finite temperature. Since the model 
exists only in Minkowski space we have worked with a for- 
mulation of the Wetterich equation in Minkowski-space 
and have Wick-rotated the momentum integrals. 

At zero temperature we have found results quite sim- 
ilar to our findings for the two-dimensional model p^ . 
An investigation of the fixed point structure yields the su- 
pcrsymmetric analogue of the Wilson-Fisher fixed-point 
for bosonic theories. It is maximally IR stable with one 
relevant direction only. As in the two-dimensional model 
the relevant direction is given by a^. Again we find a 
scaling relation between the critical exponent of 
the instable direction and the anomalous dimension rj. 
The critical exponent governs the freeze-out of the 
minimum of the superpotential W. We also find that the 
IR limit of the three-dimensional Wess-Zumino model is 
given by a conformally invariant theory. The critical ex- 
ponent governs the vanishing of the mass with the 
RG scale. This is different from three-dimensional scalar 
0{N) theories since for the supersymmetric models the 
relevant coupling does not feed back in the flow equa- 
tions of higher n-point functions. As in two dimensions 



dependently of the initial condition Aa. Nevertheless A a can be 
used to classify different theories at finite temperature. 
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we find that supersymmetry breaking is governed by the 
relevant direction. 

The finite-temperature flow equations are obtained by 
replacing the integration over the time-like direction by a 
sum over Matsubara frequencies. As bosons and fermions 
obey different thermal boundary conditions, finite tem- 
peratures introduce a soft supersymmetry breaking. 

In the phase with broken supersymmetry (restored TLi 
symmetry) at zero temperature the three-dimensional 
Wess-Zumino model flows to a masslcss field theory. In 
the ^2 symmetric phase at flnite temperature the model 
therefore behaves like a three-dimensional gas of massless 
bosons. We have found that this is indeed the case. The 
small deviations from the ideal-gas law found in our nu- 
merical studies originate from the self-interaction of the 
bosons. 

At high temperatures, supersymmetry breaking mani- 
fests itself in the fact that the flow equations for the su- 
perpotential, derived from the bosonic and the fermionic 
part of the action, are different. We also observe dimen- 
sional reduction in the way that, after a suitable rescal- 
ing, the flow equation for the "bosonic" superpotential 
in three dimensions reduces to the flow equation in two 
dimensions. Due to the absence of a fermionic thermal 
zero-mode the fermions do not contribute to the RG flow 
for small scales k <^ T. We have argued that the radius of 
convergence for the polynomial expansion of the super- 
field interpolates between the values for two and three 
dimensions as the temperature is raised. 

Even though supersymmetry is explicitly broken at fi- 
nite temperature, the Z2 symmetry of the model can be 
either restored or broken at finite temperature. Whether 
Wi2 symmetry is broken or not depends on the tempera- 
ture (and parameters of the model, i. e. the initial values 
of the couplings at the initial RG scale) . Since supersym- 
metry and Wi2 symmetry are intimately linked we have 
argued that a study of Z2 symmetry may be used to 
measure the strength of supersymmetry breaking. We 
have computed the phase diagram for the restoration of 
Z2 symmetry at finite temperatures. We find two dif- 
ferent types of phases which are separated by a second 
order phase transition: one phase with soft supersym- 
metry breaking due to finite temperature but broken Z2 
symmetry and one with restored Z2 symmetry. 

Throughout this paper we have addressed several sim- 
ilarities and differences between scalar 0{N) models and 
the Af ~ 1 Wess-Zumino model at zero and finite temper- 
atures, e. g. the fixed-point structure at zero temperature 
and the behavior at finite temperature. A detailed ex- 
ploration of both models with respect to their similarities 
and differences, in particular with respect to the large- 
limit, is deferred to future work. 
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Appendix A: Derivation of the flow equations in 
Minkowski space 

In this Addendum we derive the Wctterich fiow equa- 
tion in Minkowski space. For the sake of simplicity we 
consider a real scalar field in this appendix. The gener- 
alization to other fields, such as fermion or gauge fields, 
is straightforward. 

The generating functional in Minkowski-space is given 

by 



Z[J] ^ j Vipe 



where J denotes the external source and (J, if) = 
J d'^x J{x)lp{x). The generating functional W for the 
connected two-point functions, the so-called Schwinger 
functional, readj^ 

W[J] = i In Z[j] . 
From this we obtain 

The effective action is the Legendre transform of the 
Schwinger functional, 

r[0] = -T^[J]-(J,0), 

where is the mean (classical) field. Using ^iy[J] = 



we obtain the equation of motion for the field 



5n 



,aJW[J] SJiy) 
u y- 



SJ{y) 6(l){x) J " "(S(/)(.T 
= -Jix). 

The scale-dependent generating functional is defined as 



with 



d^q 



v{~q)Rk{q)v{q) ■ 
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The generating functional W of connected two-point functions 
should not be confused with the super potential in the main 
part. 
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Here, the momentum-dependent regulator function Rk 
provides an IR cutoff for all modes and has to sat- 
isfy three conditions: (i) Rk{p)\p^ /k^^o > which im- 
plements the IR regularization, (ii) Rkip)\k^ /p'^^o ~ 
which ensures that the regulator vanishes for fc — > 0, (iii) 
Rkip)\k^A^oo — > oo which ensures that the path integral 
is dominated by the stationary point for fc — > A — )■ oo. 
Different choices for Rk define different RG trajecto- 
ries manifesting the RG scheme dependence, but the IR 
physics Tk^o ~^ T should remain invariant provided the 
truncation captures all relevant operators for the physi- 
cal observables under investigation. In turn, a variation 
of the regulator function may lead to more insight on the 
truncation dependence of our results. 

Next, we define the scale-dependent effective action: 



-Wk[J]- J d^xJ<j)-ASk[^]. 



In order to properly formulate the flow equations in 
Minkowski-space we have to take fc^ ^ P^iP^ as flow 
parameter. Therefore we define dt = 2k^dk2 to be the 
derivative with respect to RG 'time' t — Infc^. Taking 
the derivative of Tk [0] with respect to t yields 

dtTkicI)] = -dtWkiJ] ~ dt f d^xJ(j> - dt/^Sk[<P] , (Al) 



where we have used that the source is independent of fc. 
For the derivative of Wk we obtain 



dtWk[J]^idt\nZ[J] 



Z[J] I 



(2^) 



1 



d'^q 



2Z[J] J (27r) 
1 f d'^q 



2Z[J] 



S^Z[J] 



(2n 



lidtRk)- 



6J6J 



Using the definition of the Schwinger functional, Zk [J] 



, we find 



dtWk[J]=e''^''\ J 



d'^q 



(dtRk) 



(27r)'' 5JSJ 
Now we rewrite the integrand by making use of 



(A2) 



6J5J 6J ^ ' SJ 



-iWk 



5Wk 



-iWk 



SJ ' ' 5J 
then Eq. (|A2[) can be rewritten as follows: 

1 f d'^q 



dtWk[J] 



{2ny 



-dtRk 



SWkSWk^^PwA 
5 J 5 J ^SJSJJ 



d'^q 
{2ttY 



-0 
5^Wk 



dtRk 



SJSJ 



With this relation the variation of the effective action 
Eq. (|Aip takes the form 



dtT 



(2^) 



SJSJ 



The second functional derivative of W with respect to the 
source J can be written in terms of the effective action: 



S^Tk 



SJ 



Rk 



5J ^ fS''Wk\ 
S(j) \SJSJ J 

Making use of 



Ri 



(A3) 



S{q-q') 



S SWk 



d'^q S^WkSJ 



S(l){q') 5(t> SJ J (27r)^ SJSJ S(l> ' 
we obtain the Wctterich equation in Minkowski-Space: 



dtT 



^Tr 
2 



{dtRk){ 



S'Tk 



Ri 



(A4) 



Appendix B: Derivation of the flow equations at 
flnite temperature 

In order to preserve supersymmetry in the RG flow 
for vanishing temperature, we must choose a regulator 
function which regularizes the theory in the time-like di- 
rection and the space-like directions in the same way. In 
order to make apparent how soft SUSY-breaking due to 
finite temperature emerges, we use the same regulator for 
our finite-temperature and zero-temperature studies. It 
is given by: 



r2 



\P\ 



-10 



- 1 



ri = 0. 



In the LPA, the finite-temperature flow equations can 
be obtained straightforwardly from the zero-temperature 
flow equations by replacing po by the Matsubara modes 
f„ and ujn of fcrmion and boson fields, respectively, and 
replacing the integration over pq by a summation over 
the Matsubara modes. The contribution of the bosons to 
the RG flow of our Wess-Zumino model then reads: 



jjsj' dkWi^wi!'T j -i 



d^Ps {w;!^ -p)0{e 



Ps 



where Ps denotes the momenta in space-like directions. 
Along the lines of, e.g., [4l|, we use Poisson's sum for- 
mula, 

E E / dqf{q)cM~2T:Uq), 



13 



in order to obtain 



oo 



dkWL =w:."T ^ 



oo 



dq 



— oo 



— oo 

//2 _ l.2\Qll.'2 _ J2 



87r2 

(fc2 +I^^'2)2^p2 + {2'KqTf 

For the computation of the three-dimensional integral, 
we first substitute q' — iirTq and introduce spherical 
coordinates 

p\ = r cos § sin (f, pi ^ r sin i? sin (p, q' = r cos t?. 
Performing the angular integrations yields 



dr (VP^;^/'-fc') 2rsin(f ) 

8^2^ (fc2 + ^//2)2 £ 



Finally the integration over r leads to 



dkW^ 



(M^f . ^ l-cos(f) 



87r2(fc2 + H^^'2)2 



^2 



Since we made use of Poisson's resummation formula to 
rewrite the sum over the thermal modes, we are able 
to split the flow equation into a zero-temperature and a 
finite-temperature contribution: 



dkWl 



, (W^f -fc^)W^r A,2 , .^2^1-cos(M) 

87r2(fc2 + PT/^'2)2 



1=1 



£2 



The contribution of the fermions to the RG flow of our 
model can be obtained along the lines of our derivation 
of the bosonic contribution and reads: 



87r2(fc2-f )2 



Stcrm(T) 



Introducing the dimensionless temperature T = T/k, the 
functions ghosiT) and 5fcrm(T') can be written in terms 



of Polylogarithms 
2 

/ 

2 



5bos(T) 2 [tt^ - 3Li2(e-'/^) - 3Li2(e'/^; 



5fcrm(T) = - ^fc2y2 ^ ^ 6Li2(-e-'/^) + 6Li2(-c'/^] 
Using the identity [i^ 



Li2 i'z) + Li2 (^-i^ = 2Li2(-l) - i \n{z) 
the function ghos(T) simplifies further to 

<?bos(T) = f 2 p + 1,,2 ^ _ cxp(z/f )) 

= ttT (ttT - (2si3 + if-nT + (2sb + l)2fc) - k^ , 

(Bl) 

where we have used that 

In ^exp(i/f' + infj = ^ - ztt {2sb + 1), 



SB < < sb + 1 SB ^ 



2ttT 



Similarly, exploiting the relation 



In (cxp( j/f) 



T 



2i7r sj7, 



1 1 1 

Sp — - < ^ < Sp + - => Sp = 

2 - 27rT - 2 



k 



2ttT 2 



we end up with the result 



.gform(T) ^ -k 1 - s 



27rrV 



for the fermions. As expected, the functions ghosiT) and 
5ferm(T) exhibit the same behavior as the threshold func- 
tions discussed in Ref. l25ll. 
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